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Abstract 

Analytical  expression  of  electric  power  was  deduced  in  case  of  the  large-scale  thermoelectric  device  that  consists  of  the  cylindrical 
double  tubes  like  roll  cake  and  is  exposed  to  the  two  thermal  fluids.  The  output  powers  of  16  systems  were  mathematically  described  by 
the  simultaneous  equations  based  on  heat  exchange.  The  temperature  profiles  in  the  device  depend  on  the  flow  directions  of  hot  and  cold 
fluids,  but  the  flow  directions  did  not  change  the  output  power.  Resultantly,  eight  sets  of  solutions  for  the  output  power  were  deduced.  The 
maximum  output  power  was  the  largest  in  the  two  systems  (V2CC-I  and  V2CC-II  system),  where  two  fluids  flow  in  counter  directions  and 
one  of  the  fluids  goes  into  the  system  from  the  inside  of  the  inner  cylinder.  These  chosen  systems  can  generate  the  thermoelectric  power 
equivalent  with  the  single  cylinder  system  (VIC  system)  using  only  36%  material  of  VIC. 

©  2004  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Thermoelectric  generation;  Thermoelectric  device;  Heat  transfer;  Thermal  fluid;  Bismuth  telluride 


1.  Introduction 

Thermal  energy  can  be  converted  directly  to  the  electric¬ 
ity  in  a  thermoelectric  power  plant,  where  the  thermoelectric 
motive  force  (EMF)  due  to  Seebeck  effect  is  generated  pro¬ 
portional  to  a  temperature  difference  A0  between  the  hot  and 
cold  junctions.  The  larger  A0  to  all  the  thermoelectric  panels 
is  generally  given  by  the  fluids,  which  carry  the  heat  from  the 
isolated  sources.  The  hot  fluid  warms  the  hot  junctions  from 
the  surface  of  thermoelectric  modules,  and  the  cold  fluid 
cools  the  other  junctions  in  the  same  way.  Therefore,  the  heat 
is  transmitted  to  the  cold  fluid  through  the  thermoelectric 
materials.  Basing  on  the  heat  exchange,  the  previous  stud¬ 
ies  reported  that  the  maximum  of  output  power  exists  cor¬ 
responding  to  a  certain  module  size  [1-10].  This  is  because 
the  more  serial  connections  of  thermoelements  can  generate 
the  higher  voltage,  but  because  the  output  power  decreases 
both  by  the  higher  electrical  resistance  and  by  the  smaller 
A0  due  to  the  larger  heat  exchange  through  the  modules. 

When  the  thermoelectric  panels  are  connected  three-dime- 
nsionally  [6-10],  the  dimension  of  the  thermoelectric  de¬ 
vice  can  be  shortened.  For  example,  the  flat  panels  stacked 
multiply  could  save  the  space  for  the  power  generation;  the 
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same  power  could  be  generated  in  the  four-panel-system  by 
22.5%  space  of  the  single  panel  [7,8]. 

The  form  of  cylindrical  tubes  is  more  suitable  for  fluid 
circulation.  The  multiple  tube  (T)  systems  were  reported 
previously  [6,8].  Alternatively,  the  vortical  (V)  cylindrical 
system  as  shown  in  Fig.  1  circulates  the  fluids  vortically  in 
the  gaps  between  the  thermoelectric  panels.  This  V  system 
can  flow  the  larger  amount  of  fluid  with  the  less  pressure 
loss  than  T  systems,  and  it  has  the  potential  to  produce  a 
large  power.  Fig.  1(a)  looks  like  roll  cake  and  consists  of  the 
double  thermoelectric  cylinders,  among  which  a  hot  fluid 
and  a  cold  fluid  flow  as  counter  flow. 

The  mathematical  expression  of  its  output  power  was 
previously  reported  in  the  three  models  simplified  as 
Fig.  1(b)— (d),  and  the  double  cylinder  system  with  counter 
flow  (V2CC  system,  Fig.  1(d))  was  enough  efficient  for 
miniaturization,  comparing  with  the  single  cylinder  (VI) 
systems  [10].  However,  the  other  possible  flow  directions, 
as  shown  in  Fig.  2(b)-(p),  may  hold  the  same  potential,  al¬ 
though  they  were  not  analyzed.  For  the  practical  operation 
of  the  thermoelectric  power  generation  with  the  multiple 
cylinders,  the  choice  of  fluid  directions  should  be  clarified 
from  the  fundamental  bases. 

Mode  I  and  II  in  Fig.  2  show  a  pair  of  alternative  selection 
of  the  fluid  path  whether  the  hot  fluid  flows  between  two 
thermoelectric  cylinders  or  outside  them.  Mode  I  and  II  show 
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(b)  VIC 


(a)  roll  cake  model 
(V2CC-I) 

Adiabatic  tubes 


Thermoelectric  tubes 


(c)  VIP 


(d)  V2CC-I 


Fig.  1.  Vortical  systems  with  the  cylindrical  thermoelectric  tubes. 


the  cases  that  the  fresh  fluid  enters  from  the  inside  of  the 
cylinder,  while  it  enters  from  the  outside  in  Mode  III  and 
IV.  In  Mode  I  and  III,  a  hot  fluid  is  sandwiched  by  two  cold 
fluids,  while  a  cold  fluid  is  sandwiched  by  two  hot  fluids  in 
Mode  II  and  IV. 

The  purpose  of  this  work  is  to  deduce  the  mathematical 
expressions  of  the  thermoelectric  power  extractable  from  the 
double  cylinder  systems  (Fig.  2),  and  to  search  for  the  fur¬ 
ther  improvements  in  the  operation.  The  choices  of  flow  di¬ 
rections  and  fluid  paths  suitable  for  high  output  are  focused. 


2.  Basic  assumptions  and  models 

A  brief  outline  of  the  methods  for  modeling  and  mathe¬ 
matical  analysis  will  be  given  here  because  their  details  of 
the  common  bases  were  reported  separately  [6-10]. 

2.7.  Directions  and  paths  of  fluids 


and  three  paths  for  fluids.  The  two  fluids  occupy  these  paths, 
and  one  of  the  fluids  is  circulated  twice.  The  systems  are 
classified  by  the  direction  of  fluid  flow  (counter  flow,  “C”, 
and  parallel  flow,  “P”)  and  the  choice  of  fluid  paths  (Mode 
I-IV).  For  example,  in  the  V2CC-I  system,  the  cold  fluid 
flows  the  most  inner  path  (Mode  I)  and  all  the  fluids  in  three 
paths  flow  in  counter  mode  (first  C  for  inner  cylinder  and  sec¬ 
ond  C  for  the  outer  cylinder).  The  hydrodynamic  boundary 
conditions  for  1 8  mathematical  models  are  listed  in  Table  1 . 

2.2.  Cylindrical  thermoelectric  tubes 

Our  thermoelectric  module  consists  of  cylindrical 
fan- shaped  thermoelements  with  a  single  layer  of  II -type 
p-n  junctions,  as  shown  in  Fig.  3  [10].  All  the  thermoele¬ 
ments  are  homogeneously  aligned  and  connected  electri¬ 
cally  in  series.  The  hot  and  cold  fluids  flow  along  the  both 
cylindrical  surfaces  in  the  direction  of  ±cp  (see  Fig.  1(d)). 
The  temperature  of  the  fluid  j,  Tj(cp),  becomes  as  a  func¬ 
tion  of  <p .  The  inhomogeneities  such  as  fins,  electrodes  and 
isolators  are  neglected  or  considered  as  the  heat  transfer 
coefficient,  h ,  of  the  smooth  surface. 


3.  Equations  for  output  power 

3.1.  Deduction  of  equations 


Considering  the  angular  dependency  of  T )(sp),  the  heat 
transfer  coefficient  K[  is  first  calculated  as  the  over-all  heat 
transfer  coefficient  per  unit  length  through  the  cylinder  i  in 
the  direction  perpendicular  to  r  axis  [9-12]: 

1 

Ki  =  - .  (1) 

(1/ h i, i i )  +  lfr(n,i+l/n,i)/^  T  1/ 


where  rz- j  and  hij  are  the  radius  of  the  cylinder  and  the  heat 
transfer  coefficient  between  the  cylinder  i  and  the  fluid  j ,  re¬ 
spectively.  A  is  the  average  heat  conductivity  of  the  module, 
defined  by 


A  = 


^pTp  +  AnLn  ApTp  +  AnLn 


+  Ln 


L 


where  Ap,  An,  Tp  and  Ln  are  the  heat  conductivity  and  the 
length  of  p-  and  n-type  elements,  respectively.  Fig.  3  shows 
the  assumed  setup  of  thermoelectric  junctions.  L  is  the  unit 
size  of  a  thermoelectric  pair.  Therefore,  L  =  Lp  +  Ln. 

Because  our  systems  consist  of  two  cylinders,  a  set  of 
three  simultaneous  derivative  equations  is  written  from  the 
heat  transfers  through  two  cylinders: 

d7j  (cp) 

MiCp,i— ^  =  ±KlW(Ti(<p)  -  T2(cp )),  (3) 

acp 


For  easy  mathematical  handling,  the  roll  cake  systems 
(Fig.  1(a))  are  simplified  to  the  systems  with  the  combined 
multi-cylinders  (Fig.  2).  Therefore,  the  vortical  double  sys¬ 
tems  (V2  systems)  consist  of  two  thermoelectric  cylinders 


d  T2(<p) 

±M2Cp, 2— =  ±Klw(Tl((P)  -  T2(cp)) 
acp 

±K2w(T2(cp)  -  Tflcp)),  (4) 
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(m)  V2CC-IV 


(n)  V2CP-IV 


(o)  V2PP-IV 


(p)  V2PC-IV 


(a)  V2CC-I 


(b)  V2CP-I 


(c)  V2PP-I 


(d)  V2PC-I 


(e)  V2CC-II 


(f)  V2CP-II 


(g)  V2PP-II 


(h)  V2PC-II 


(i)  V2CC-III 


(j)  V2CP-I1I 


(k)  V2PP-II1 


(1)  V2PC-III 


Fig.  2.  Analyzed  systems  with  the  double  cylindrical  thermoelectric  tubes. 


d  73  0) 

±M3CP,3— =  ±K2w(T2(cp)  -  T3(<p)),  (5) 

d  (p 

where  Mi  and  Cpj  are  mass  flow  rate  and  heat  capacity 
for  fluid  i,  respectively.  For  simplicity,  the  steady  state  is 


i 


Fig.  3.  Dimensions  of  thermoelectric  modules. 


considered,  and  the  heat  transfer  by  Peltier  effect  and  the 
exothermal  heat  by  Joule  effect  are  neglected.  The  choice 
of  sign  ‘zb’  in  Eqs.  (3)— (5)  depends  on  the  system  design. 
The  temperature  Oij(<p)  at  the  p-n  junction  in  the  cylinder  i 
facing  to  the  fluid  j  is  written  using  Tj(<p): 


0i,i+\{<p)  -  Oi,i(<p) 


A 

X 


[Ti+xiv) 


Ti(fp)}  In 


D',z'+ 1 


ri,  i 


The  electromotive  force,  E ,  is  the  summation  of  EMF  over 
the  two  cylinders: 


E  —  E\  +  E2 


N(f)nxwa 

2n 


Oi+i,i(<p)}d<p  , 


where  a  is  the  difference  between  the  Seebeck  coefficients 
of  p  and  n  elements,  and  nx  are  the  number  of  thermo¬ 
electric  pairs  along  a  cylindrical  circulation,  and  the  number 
density  of  pairs  in  the  cylinder  length  direction,  respectively. 
Therefore,  N $  =  27r/</>  and  nx  =  1/L.  0  is  the  circumfer¬ 
ential  angle  of  a  pair.  The  electric  resistance,  Ri,  is  deduced 
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Table  1 


Conditions  for  the  analyzed  systems 


System 

Direction  of  fluidsa 

Boundary  conditions 

Model 

Reference 

VIC 

d\  >  0,  d2  <  0 

Ti(0)  =  7*\  T2(2tt)  =  Tf 

Fig.  1(b) 

[10] 

VIP 

d\  >  0,  d2  >  0 

Ti( 0)  =  Tf,  UCO)  =  Tf 

Fig.  1(c) 

[10] 

V2CC-I 

d\  >  0,  d2  <  0,  d3 

—  d\  >0 

Ti(0)  =  Tf,  72(2tt)  =  Tf,  73(0)  =  7i(2tt) 

Figs.  1(d)  and  2(a) 

[10] 

V2CC-II 

ibid. 

Ti(0)  =  Tf,  T2(2tc)  =  Tf,  73 (0)  =  7i(2tt) 

Fig.  2(e) 

V2CC-III 

ibid. 

Ti(0)  =  73(2tt),  T2(2tc)  =  Tf,  T3( 0)  =  Tf 

Fig.  2(i) 

V2CC-IV 

ibid. 

Ti(0)  =  T3(2tc),  T2(2jc)  =  Tf,  T3( 0)  =  Tf 

Fig.  2(m) 

V2CP-I 

d\  >  0,  d2  <  0,  d3 

=  —d]  <  0 

T\(0)  =  Tf,  T2(2jt)  =  Tf,  T3(2tt)  =  Tx(2n) 

Fig.  2(b) 

V2CP-II 

ibid. 

Ti( 0)  =  Tf,  T2(2jt)  =  Tf,  T3(2tt)  =  Tx(2n) 

Fig.  2(f) 

V2CP-III 

ibid. 

Ti(0)  =  r3 (0),  T2(2tt)  =  Tf,  T3(2tt)  =  Tf 

Fig.  2(j) 

V2CP-IV 

ibid. 

Ti(0)  =  7M0),  T2(2tt)  =  Tf,  T3(2tt)  =  Tf 

Fig.  2(n) 

V2PP-I 

d\  >  0,  d2  >  0,  d3 

—  d\  >  0 

Ti(0)  =  Tf,  7M0)  =  Tf,  73(0)  =  T\(2n) 

Fig.  2(c) 

V2PP-II 

ibid. 

Ti(0)  =  Tf,  72(0)  =  Tf,  73(0)  =  T\(2j[) 

Fig.  2(g) 

V2PP-III 

ibid. 

71(0)  =  T3(2tc),  72(0)  =  Tf,  73(0)  =  Tf 

Fig.  2(k) 

V2PP-IV 

ibid. 

Ti(0)  =  T3(2jc),  72(0)  =  Tf,  7M0)  =  7jn 

Fig.  2(o) 

V2PC-I 

d\  >  0,  d2  >  0,  d 3 

=  —d\  <  0 

Ti(0)  =  Tf,  72(0)  =  7^n,  T3(2tt)  =  Tx(2n) 

Fig.  2(d) 

V2PC-II 

ibid. 

Ti( 0)  =  Tf,  72(0)  =  Tf,  T3(2jt)  =  T^n) 

Fig.  2(h) 

V2PC-III 

ibid. 

Ti(0)  =  T3( 0),  7V0)  =  Tf,  T3(2tt)  =  Tf 

Fig.  2(1) 

V2PC-IV 

ibid. 

Ti(0)  =  r3(0),  7V0)  =  Tf,  T3(2tt)  =  T™ 

Fig.  2(p) 

a  When  the  vector  of  fluid  i  is  defined  as  M;  at  (p  =  0,  M;  is  written  as  M/  =  diet,  where  et  is  an  unit  vector  in  the  positive  direction  at  <p  =  0  and 
di  is  a  constant.  The  mass  flow  rate,  Af/,  is  equal  to  I  <7/1. 


by  integrating  the  resistance  of  a  thin  fan-shaped  layer  in 
the  direction  of  radius: 


=  /  Pp  Pn\  (ru+l\ 

'  2n  Up  L»)  n  V  n,  ) 


where  pv  and  pn  are  the  specific  resistivity  of  p  and  n  ele¬ 
ments,  respectively. 

Finally,  the  output  power,  P ,  is  optimized  by  balancing 
the  internal  and  external  resistances: 

E 2 

P= - .  (9) 

MRi  +R2) 


This  work  uses  hereafter  this  P ,  where  the  internal  resistance 
is  equal  to  the  external  resistance. 


3.2.  Equations  of  output  power 


The  simultaneous  equations  Eqs.  (3)— (5)  could  be  ana¬ 
lytically  solved,  but  the  explicit  expression  of  the  solutions, 
Tj((p),  were  complicated  and  the  mathematical  form  of  P  be¬ 
came  very  long.  The  following  conditions  are  used  hereafter 
in  order  to  express  simply  the  analytical  difference  among 
the  V2  systems: 


M\Cpp\  =  \M2Cpp\  =  \M2Cp^\  =MCp. 


(10) 


The  effect  of  MCP  on  P  was  discussed  for  V 1  systems  [10]. 

Under  these  situations,  the  simultaneous  equations 
Eqs.  (3)— (5)  can  be  re-written  as 


d  /  Ti(<p)\  /  Ti(cp)\ 

—  T2(<p)  \=m\  T2(cp)  .  (11) 

d(p\T3(<p))  \T3(<p)J 

where  the  matrix,  m,  depends  on  the  system  design  as  listed 
in  Table  2. 

Using  the  procedure  mentioned  above,  the  output  power 
for  the  V2  systems  was  calculated  and  listed  also  in  Table  2. 
The  definitions  of  Ka i,  Ka2  and  X  are  identical  with  the 
previous  works  [9,10].  The  result  for  V2CC-I  system  agreed 
with  the  previous  report  [10].  The  temperature  dependencies 
of  all  the  variables  are  neglected  for  simplification.  Although 
they  can  affect  the  output  power,  the  output  power  cannot  be 
deduced  in  the  analytical  way  if  we  consider  the  temperature 
dependencies  of  the  materials’  parameters.  The  numerical 
solution  with  such  considerations  will  be  reported  separately. 

First  it  is  noted  that  the  output  powers  of  Mode  I  were 
equal  to  those  of  Mode  II  in  four  types  of  fluid  direction. 
Similarly  those  of  Mode  III  were  equal  to  those  of  Mode  IV. 
Namely,  the  path  choice  for  the  hot  fluid  did  not  affect  the 
output  power.  Only  a  difference  was  the  sign  of  EMF. 

Second,  the  output  power  contained  the  common  term  of 
A  =  nxM2Cp2LnLpa2(fn  -  T/)2 /kw(LpPn  +  LnPp).  The 
output  power  was  simply  proportional  to  nxa2(TE'  —  T™)2 . 
All  the  output  powers  were  independent  from  the  number 
of  thermoelectric  elements,  Nq.  The  contributions  of  the 
other  variables  are  not  simple  because  they  affect  through 
the  other  terms. 

The  other  parts  were  the  complicate  functions  of  sys¬ 
tem  geometry,  material  properties  and  overall  heat  transfer 


Table  2 

Matrix  for  the  differential  equation  and  the  output  power,  where  M\CPfi  =  \M2CPt2 |  =  \M2CP^\  =  MCP 


System 


Matrix,  m 


Output  power,  P 


V2CC-I  and  V2CC-II  m  = 

V2CC-III  and  V2CC-IV  ibid. 

V2CP-I  and  V2CP-II  m  = 

V2CP-III  and  V2CP-IV  ibid. 

V2PP-I  and  V2PP-II  m  = 

V2PP-III  and  V2PP-IV  ibid. 

V2PC-I  and  V2PC-II  m  = 

V2PC-III  and  V2PC-IV  ibid. 


(t 


(v 


(I 


(I 


fi  0  \ 

h  +  h  -h  p 
h  -h ) 


p 


fi  o  \ 
fi  +  h  -h  p 
-h  h 


p 


h  o  \ 
-/i  -h  h  I  P 
h  -h ) 


p 


h  0  \ 

-/I  -h  h\  P 
-h  /2  / 


p 


A 


2(Ka\  +  Ka2 ) 


'sinh(X){(^i^ii  +  K2Ka2)  cosh(X)  +  jK^K2{Kax  -  Ka2 )  sinh(X)} 
JWkT2  cosh(2X)  +  (Ki  +  K2)  sinh(2X) 


-.2 


A 


2(Ka\  +  Ka2 )  |_ 
A 

%{Ka\  +Ka2 )  |_ 


'sinh(X){(^ifiai  +  K2Ka2 )  cosh(X)  +  JKfKi{-Kai  +  Ka2 )  sinh(X)} 
jKfT2  cosh(2X)  +  (Ki  +  K2)  sinh(2X) 


-.2 


A 

8C Ka\  +Ka2 )  |_ 
A 

8(X<*i  +  Xa2)  L 

A 

8(X<*i  +  Xa2) 

A 

8(X<zi  +Ka2)  |_ 
A 

+  Xa2)  L 


^K^(Kai  -  Ka2){cosh(2Z)  -  F}  +  V^l  +  X2Xa,  sinh(2Z) 
yZ^cosh(2Z)  +  V^i  +  Kl  sinh(2Z) 

Vk>(Kai  -  Ka2)(cosh(2Z)  -  l/F)  -  *JKX  +  K2Kci\  sinh(2Z) 

cosh(2Z)  +  y/Ki  +  K2  sinh(2Z) 


n  2 


-.2 


n2 


MCpU{Ka\  exp(— 2V)  —  Ka2  exp(2V)  —  (. Ka\  —  Ka2 )  cosh(2£/)}  +  nw(2K\Ka\  —  K2Ka\  —  K\Ka2  +  2K2Ka2)  sinh(2 U) 

2 U  exp(2V)  +  U  cosh(2 U)  +  V  sinh(2 U) 

MC pU {—Ka\  exp(— 2V)  +  Ka2  exp(2V)  +  (. Ka\  —  Ka2 )  cosh(2t/)}  +  nw(2K\Ka\  —  K2Ka\  —  K\Ka2  +  2K2Ka2)  sinh(2£/) 

2 U  exp(2V)  +  U  cosh(2 U)  +  V  sinh(2 U) 

\fK[{Kai  -  Ka2){cosh(2Y)  -  E]  +  y/Ki  +  K2Ka2  swh(2Y)~'2 

cosh(27)  +  y/Ki  +  K2  sinh(27) 


n  2 


yfK[{Ka\  -  X<32){cosh  (27)  -  1/E}  -  y/Ki  +  K2Ka2  sinh(2 7) 
_ *JK\  cosh(27)  +  JKi  +  K2  sinh(2 7) _ 


-1 2 


Where  f\  =  K\w/MCp,  f2  =  K2w/MCp,  A  =  nxM2Cp2LnLpa2(T[n  -  T1^)2 /wk(Lppn  +  Lnpp),  E  =  exp(-2^/i),  F  =  exp(27i/2),  £/  =  W/i2  “  /1/2  +  /22,  V  =  7r(/i  +  /2),  X  =  njjljji, 
y  =  ttVTTCA  +  /2),  z  =  ^ V./2 (./1  +  /2b  =  (1/ 7tA)  ln(ri2/rn),  and  Ka2  =  (l/nX)  ln(r23/r22). 


to 
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coefficient  K[.  Because  Kf  depended  on  the  radius  as  written 
in  Eqs.  (1)  and  (2),  our  expressions  could  not  be  simplified 
analytically  by  the  similar  way  with  the  flat  multi-panels 
systems  [7,8].  Therefore,  the  output  power  will  be  numeri¬ 
cally  compared  using  the  thermo-physical  values  of  Bi2Te3 
semiconductors. 


4.  Numerical  comparison  of  systems 

4.1.  Physical  properties  and  module  optimization 


Table  3  shows  a  typical  set  of  properties  for  an  identical 
sample  of  the  BiTe  thermoelectric  element  at  room  temper¬ 
ature  [13],  and  used  for  this  work.  Table  4  shows  the  fluid 
properties  and  the  parameters  for  thermoelectric  cylinders, 
considering  the  utilization  of  nitrogen  gas.  The  heat  transfer 
coefficient,  h ,  was  set  to  be  independent  of  the  flow  rate  and 
of  the  cross-section  of  fluid  path.  These  values  are  identical 
with  our  previous  studies  [6-10]. 

The  output  power  P  in  Table  2  was  optimized  against  the 
geometrical  factors,  Ln  and  Lp,  using  that 


Ln  = 


-y/^pAl 


+  yApPn 


L  and  Lp  = 


^/An/9p 


y/^nPp  T  y^A-pPn 


L 


(12) 


Eq.  (12)  is  generally  taken  for  the  analysis  with  the  infi¬ 
nite  heat  sources,  and  it  was  also  a  good  approximation  for 
cylindrical  geometry  [9,10]. 


4.2.  Maximum  output  power 


Fig.  4  shows  18  kinds  of  the  output  power,  which  were 
thus  optimized  by  Eq.  (12)  and  calculated  by  the  data  given 


Fig.  4.  Comparison  of  the  output  power  of  the  systems,  where  the  cylinder 
length  w  =  2  m. 


in  Tables  3  and  4.  Two  solutions  for  V 1  system  are  referred 
from  [10].  The  output  powers  of  Mode  I  and  III  exactly 
agreed  with  those  of  Mode  II  and  IV,  respectively,  naturally 
because  the  mathematical  expressions  are  identical  as  shown 
in  Table  2. 

As  the  inner  radius,  rn,  increased,  the  calculated  output 
power  increased  to  a  maximum,  Pmax,  at  an  optimized  radius, 
r^1',  and  then  decreased  monotonically.  This  phenomenon 

was  common  for  all  systems,  but  Pmax  and  were  charac¬ 
teristic  for  the  individual  systems.  It  was  difficult  to  deduce 


Table  3 

Specific  values  of  thermoelectric  elements  used  here 


Materials 

Seebeck  coefficient,  a  (pV7K) 

Resistivity,  p  ([x£2m) 

Thermal  conductivity,  X  (W/Km)  Figure  of  merit,  Z  (K  1 ) 

Bi-54.3  at.%  Te  (p) 
Bi-64.5  at.%  Te  (n) 

162 

-240 

5.55 

10.1 

2.06  2.605  x  10"3 

2.02 

Table  4 

Parameters  for  thermoelectric  power  generation  system 

Variables 

Values  used  in  this  work 

Thermoelectric  device  Length,  w 

Thickness  of  tube,  d 
Inner  radius  of  inner  tube,  r\\ 

Distance  between  two  tubes,  a  (fluid  path  height) 
Number  density  of  pairs,  nx 

Number  of  pairs  round  a  circumferential  cycle,  N. </> 


Variable  (typically  2  m) 
Variable  (typically  0.05  m) 
Variable 

Variable  (typically  0. 1  m) 
50m_1 
100 


Thermal  sources 


Hot  source 
Cold  source 


N2  gas  (inlet  T jn  =  500  K) 
N2  gas  (inlet  Tcin  =  300  K) 


Gas  properties 


Specific  heat,  Cp  \i Qt  —  Cp,coo l 
Mass  flow  rate,  Mh0t  =  Mcooi 

Heat  transfer  coefficient,  h\  \  =  h\  2  =  ^2,2  =  ^2,3  =  ^hot  =  hcoo\ 


1044  J/kgK  (at  400  K,  0.1  Pa) 

1.00  kg/s 

500  W/m2  K  independent  of  flow  rate 
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Table  5 

Maximum  output  power,  Pmax ,  and  thermoelectric  motive  force,  EopL ,  at  the  optimized  radius,  r°fi^ 


System 

Length  of  the  tubes,  w 

1  m 

2m 

5  m 

pmax 

£°pt.  (V) 

‘  (m) 

pmax 

£°pt-  (V) 

'  (m) 

pmax  ^yy) 

£°pt-  (V) 

r°u'  (m) 

V2CC-I  and  V2CC-II 

5531 

375.6 

1.880 

5594 

737 .4 

0.903 

5703 

1812 

0.3345 

V2CC-III  and  V2CC-IV 

5361 

362.2 

1.883 

5256 

714.0 

0.905 

4881 

1670 

0.3375 

V2CP-I  and  V2CP-II 

5231 

379.8 

1.661 

5273 

759.3 

0.793 

5302 

1845 

0.2922 

V2CP-  III  and  V2CP-IV 

5070 

373.8 

1.662 

4954 

735.1 

0.795 

4539 

1702 

0.2946 

V2PP-I  and  V2PP-II 

5010 

389.5 

1.504 

5056 

778.4 

0.716 

5084 

1890 

0.2611 

V2PP-  III  and  V2PP-IV 

4855 

383.3 

1.505 

4749 

753.9 

0.717 

4349 

1744 

0.2627 

V2PC-I  and  V2PC-II 

5254 

380.7 

1.660 

5320 

762.7 

0.793 

5424 

1870 

0.2908 

V2PC-  III  and  V2PC-IV 

5092 

374.6 

1.662 

4997 

738.3 

0.795 

4636 

1725 

0.2925 

pmax  an(j  ropt-  anaiytically,  because  they  depended  on  heat 
transfer  coefficient  K[  in  the  complicated  form.  Pmax  and 
were  numerically  solved  using  the  equations  in  Table  2, 
and  they  are  listed  in  Table  5. 

4.3.  Temperature  profiles 

Figs.  5  and  6  show  the  temperature  distributions  when  the 
cylinder  radius  was  equal  to  the  optimum  inner  radius,  r° 
The  surface  temperatures,  Oij((p ),  are  shown  in  addition  of 
the  fluid  temperatures,  Tficp).  The  hot  and  cold  fluid  are 
cooled  and  warmed,  respectively,  as  they  proceed  along  the 
module  surface.  The  temperature  profiles  varied  depending 
on  the  system  design. 

As  shown  in  V2CC  systems  (Fig.  5),  the  temperature 
profiles  for  Mode  I  and  II  and  for  Mode  III  and  IV  were 


symmetric,  and  this  causes  the  agreement  of  the  output 
power  in  these  modes.  Therefore,  the  temperature  profiles 
for  the  other  systems  were  shown  only  for  Mode  I  in 
Fig.  6. 

As  reported  previously  [10],  the  fluid  temperatures  of  VIC 
system  (Fig.  6)  changed  linearly  when  M\Cpp  =  \MiCpp\, 
and  the  outlet  temperatures  of  the  hot  and  cold  fluids  were 
the  same,  400  K  (7£"‘  =  =  (7^  +  7^)/ 2).  However, 

the  fluid  temperatures  of  V2  systems  were  not  linear.  The  hot 
fluid  sandwiched  by  two  cold  fluids  was  cooled  rapidly,  and 
the  temperature  rises  of  coolants  were  slow.  At  the  optimum 
conditions  for  Pmax,  the  outlet  temperatures  of  the  hot  and 
cold  fluids  (T^  and  T^d)  were  above  and  below  400  K, 
respectively.  It  means  that  the  thermal  energy  of  the  fluids 
in  V2  systems  did  not  exchanged  so  completely  through  the 
thermoelectric  panels  as  that  in  VI  systems. 


Angle,  cp  /rad 


Angle,  (p  /rad 


Angle,  <p  /rad 


Angle,  cp  /  rad 


Fig.  5.  Temperature  profiles  of  V2CC  system  along  the  cylinder  surfaces, 
when  the  output  power  becomes  the  maximum.  The  arrows  show  the  fluid 
directions.  T\,  T2  and  T$  are  temperatures  of  the  fluids,  and  0\\,  612 ,  O22 
and  $23  are  the  surface  temperatures  ( w  =  2  m). 


Fig.  6.  Temperature  profiles  along  the  cylinder  surfaces,  when  the  output 
power  becomes  the  maximum.  The  arrows  show  the  fluid  directions.  T\, 
T2  and  73  are  temperatures  of  the  fluids,  and  6\\,  612 ,  $22  and  622  are 
the  surface  temperatures  ( w  =  2  m). 
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4.4.  Largest  Pmax 

The  existence  of  Pmax  was  reported  also  in  the  flat  panels 
and  the  tubes  [1-10].  It  is  partially  because  the  larger  ther¬ 
moelectric  panel  size  reduces  the  temperature  difference  AT 
between  the  two  fluids. 

As  shown  in  Table  5  and  Fig.  4,  the  maximum  output 
power  of  V2CC-I  and  V2CC-II  systems,  Pmax(V2CC-I 
and  V2CC-II),  was  the  largest,  and  Pmax(V2PP-III  and 
V2CC-IV)  was  the  smallest  in  the  V2  systems  when  the 
cylinder  length  w  was  fixed.  However,  r^t'(V2CC-I  and 

V2CC-II)  was  the  largest  and  r°f  (V2PP-III  and  V2CC-IV) 
was  the  smallest.  This  tendency  consists  with  that  for 
VI  systems  [10].  The  reason  why  Pmax (counter  flow)  > 
pmax  (parapel  flow)  can  be  qualitatively  understood  by  the 
temperature  profiles  (Fig.  6)  that  J^ATXcounter  flow)  is 
larger  than  A  ^(parallel  flow)  [7,11,12].  AT(V2CC)  is 
nearly  constant  for  all  the  thermoelectric  elements,  while 
AT(parallel  flow)  is  large  at  the  inlet  but  it  decreases  sig¬ 
nificantly.  Therefore,  the  counter  flow  can  generate  larger 
power  than  the  parallel  flow,  although  the  system  size  of 
counter  flow  to  obtain  the  maximum  output  is  larger  than 
that  of  parallel  flow  [10]. 

It  is  noteworthy  that  Pmax(V2CC-I  and  V2CC-II)  was 
nearly  equal  to  pmax  in  VI  systems,  i.e.,  98.7%  of 
pmax(ViC)  and  117.5%  of  Pmax(VlP),  when  w  =  2m. 
However,  r^t'(V2CC-I  and  V2CC-II)  was  only  45.5%  of 

r^L(VlC)  and  61.6%  of  r^"(V IP).  The  amount  of  thermo¬ 
electric  material  necessary  for  Pmax(V2CC-I  and  V2CC-II) 
is  evaluated  as  36.0  and  91.1%  of  those  for  Pmax(VlC) 
and  Pmax(VlP),  respectively.  V2CC  system  is  thus  more 
compact  and  material  saving  than  V 1  systems. 

The  power  difference  between  Pmax(V2CC)  and 
pmax(y2pp)  was  smaller  than  that  between  Pmax(VlC)  and 
Pmax(\lP),  and  the  radius  difference  between  r^L(V2CC) 

and  r^l'(V2PP)  was  closer  than  that  between  r^t-(VlC) 

and  r^L(VlP).  This  means  that  the  choice  of  flow  direc¬ 
tions  in  V2  systems  was  not  so  sensitive  as  in  VI  systems 
to  obtain  the  large  output  power. 

4.5.  V2CP  and  V2PC  systems 

The  output  power  of  V2CP-I  and  V2CP-II  systems, 
P(V2CP-I  and  V2CP-II),  was  very  close  to  P(V2PC-I  and 
V2PC-II)  in  the  wide  range  of  radius,  rn,  as  shown  in  Fig. 
4.  Similarly  P(V2CP-III  and  V2CP-IV),  was  very  close  to 
P(V2PC-III  and  V2PC-IV).  The  combination  of  counter 
and  parallel  flow  may  cancel  their  characters  each  other, 
and  the  output  powers  of  V2CP  and  V2PC  systems  may 
become  similar. 

As  illustrated  in  Fig.  2,  these  two  systems  need  the  sudden 
change  of  flow  direction  at  <p  =  0  or  27r.  This  causes  the 
large  pressure  loss  for  fluid  circulation.  When  we  consider 
the  realistic  application  in  the  V2  systems,  therefore,  the 
vortical  continuous  fluid  flow  such  as  Fig.  1(a)  would  be 


Fig.  7.  Temperature  difference  between  the  cylinder  surfaces  of  V2CC 
systems,  when  the  output  power  becomes  the  maximum. 

completed  not  by  V2CP  or  V2PC  systems,  but  by  V2CC  or 
V2PP  systems. 

4.6.  Mode  I  and  III 

The  output  powers  of  Mode  I  and  II  were  always  larger 
than  those  of  III  and  IV,  as  shown  in  Fig.  4.  This  reason 
can  be  found  in  the  length  of  fluid  path  as  discussed  below. 
Here,  V2CC-I  and  V2CC-III  systems  are  picked  up  as  the 
example. 

The  temperature  drop  of  the  hot  fluid  was  nearly  the  same 
both  in  two  systems  as  shown  in  Fig.  5.  However,  the  tem¬ 
perature  rises  of  the  coolant  were  different.  The  surface  tem¬ 
perature  differences  AO  were  shown  in  Fig.  7.  AO  of  V2CC-I 
system  were  always  larger  than  AO  of  V2CC-III.  It  is  ex¬ 
plained  by  the  difference  of  fluid  path  length,  2nr ,  at  the 
outer  and  inner  paths.  Because  the  cold  fluid  at  V2CC-III 
system  is  warmed  during  the  first  circulation  in  the  longer 
outer  path,  73(V2CC-III,  cp  =  2tt)  becomes  higher  than 
74(V2CC-I,  (p  =  27t).  The  fluid  after  the  first  circulation 
with  this  higher  temperature  of  V2CC-III  system  enters  into 
the  inner  path  for  the  second  circulation.  Then  the  AO  of  the 
inner  cylinder  becomes  smaller.  As  deduced  in  Eq.  (7),  the 
summation  of  AO  is  proportional  to  E.  Therefore,  V2CC-I 
system  can  generate  larger  E  than  V2CC-III  system.  Note 
that  r^t'(V2CC-I)  is  nearly  equal  to  r^t'(V2CC-III),  and  the 
internal  electric  resistance  of  power  plant  is  the  same.  We 
may  conclude,  therefore,  that  the  larger  temperature  differ¬ 
ence  of  the  fluids  should  be  given  at  the  surface  of  the  inner 
thermoelectric  cylinder  for  higher  P. 

4. 7.  Effects  of  panel  thickness  and  path  width 

The  thicker  thermoelectric  cylinder  can  keep  the  hot 
fluid  temperature  high,  because  the  heat  exchange  through 
the  cylinder  delays.  Therefore,  the  higher  output  power 
is  expected  when  the  thermoelectric  cylinder  becomes 
thicker.  Using  the  analytical  expression  of  the  output  power 
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Thickness,  dim 


Fig.  8.  Optimum  thickness  for  the  maximum  output  power  in  V2CC-I 
system,  where  the  cylinder  length  w  =  2  m.  r\ i  and  a  are  the  inner  radius 
of  the  inner  cylinder  and  the  fluid  path  height,  respectively. 

from  V2CC-I  system  (Table  2)  as  a  model,  the  effect  of 
the  panel  thickness,  d ,  on  P(V2CC-I)  was  calculated  in 
Fig.  8. 

When  the  inner  radius  of  the  cylinder,  r\ i ,  was  fixed,  P  has 
the  maximum  at  an  optimum  thickness.  The  decrease  of  P 
at  the  thicker  panel  was  due  to  the  larger  electric  resistance. 
This  maximum  of  P ,  Pmax,  became  larger  as  r\\  increased. 
Simultaneously  the  optimum  thickness  for  Pmax  increased. 
This  is  due  to  the  smaller  drop  of  temperature  difference 
between  two  fluids  when  the  heat  exchange  becomes  the 
less. 

The  path  height  of  the  fluid,  or  the  gap  distance  between 
the  cylinders,  a ,  did  not  change  P  seriously,  a  and  d  did  not 
affect  the  optimum  radius,  .  However,  it  is  noted  that  this 
work  used  a  constant  heat  transfer  coefficient,  h.  Practically, 
h  depends  on  the  mass  flow  rate  and  the  cross-section  of  the 
path  [11,12].  Therefore,  the  influence  of  a  will  be  studied 
further. 


5.  Conclusion 

This  work  studied  the  thermoelectric  power  generation 
with  cylindrical  double  tubes  in  which  thermoelectric  el¬ 
ements  were  embedded.  The  fluid  path  connections  were 
designed  as  the  roll-cake.  Sixteen  kinds  of  the  system  de¬ 


sign  were  analyzed  using  the  steady  heat  balance,  and  their 
output  powers  P  were  expressed  in  eight  mathematics  equa¬ 
tions.  The  V2CC-I  and  V2CC-II  system  were  chosen  by 
the  largest  output  power  Pmax  at  the  optimum  radius.  Their 
two  cylinders  were  exposed  to  the  two  counter  flows.  Be¬ 
cause  the  temperature  difference  of  two  fluids  can  be  kept 
larger  after  the  first  circulation,  the  fresh  fluids  should  enter 
from  the  inside  of  the  inner  cylinder  (Mode  I  and  II).  V2CP 
and  V2PC  systems  may  cause  the  larger  pressure  loss,  and 
V2PP  systems  did  not  show  the  large  Pmax.  V2CC-I  and 
V2CC-II  system  can  generate  the  equivalent  power  to  VIC 
systems,  using  only  36%  material  of  VIC. 
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